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Abstract 

Let 77 = (f](t))t£T be a sample continuous max-infinitely random field on a locally 
compact metric space T. For a closed subset S £ T, we note 775 the restriction of 
77 to S. We consider f3(Si,S2) the absolute regularity coefficient between rjs 1 and 
rjs 2 , where Si, 62 are two disjoint closed subsets of T. Our main result is a simple 
upper bound for /3(Si, S2) involving the exponent measure /x of r\: we prove that 
/3(<Si> £2) < 2 /P[r/ /, 77 ^5 2 f][J>(df), where / ^5 5 means that there exists 
s £ S such that /(s) > g(s). 

If 77 is a simple max-stable random field, the upper bound is related to the 
so-called extremal coefficients: for countable disjoint sets Si and S2, we obtain 
/3(5i,5 2 ) < 4^( sl)S2 ) e s lX s 2 (2 -9(si,s 2 )), where 0(si,s 2 ) is the pair extremal coef- 
ficient. 

As an application, we show that these new estimates entail a central limit theorem 
for stationary max-infinitely divisible random fields on Z d . In the stationary max- 
stable case, we derive the asymptotic normality of three simple estimators of the 
pair extremal coefficient. 
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1 Introduction 



Max-stable random fields turn out to be fundamental models for spatial extremes since 
they arise as the the limit of rescaled maxima. More precisely, consider the component- 
wise maxima 

Vn(t) = max t G T, 

l<i<n 

of independent copies i > 1, of a random field £ = (£(i)) te 7\ If the random field 
Vn = {Vn{t))teT converges in distribution, as n — > oo, under suitable affine normalization, 
then its limit 7] = {i](t)} t £T is necessarily max-stable. Therefore, max-stable random 
fields play a central role in extreme value theory, just like Gaussian random fields do in 
the classical statistical theory based on the central limit Theorem. 

Max-stable processes have been studied extensively in the last decades and many of 
their properties are well-understood. For example, the structure of their finite dimensional 
distributions is well known and insightful Poisson point process or spectral representations 
are available. Also the theory has been extended to max-infinitely divisible (max-i.d.) 
processes. See for example the seminal works by Resnick [20], de Haan [TT1 18], de Haan 
and Pickands |10] . Gine Hahn and Vatan [16], Resnick and Roy [21] and many others. 
More details and further references can be found in the manographs by Resnick [20J or 
de Haan and Fereira [9]. 

The questions of mixing and ergodicity of max-stable random process indexed by M 
or Z have been addressed recently. First results by Weintraub [28] in the max-stable 
case have been completed by Stoev [25] . providing necessary and sufficient conditions 
for mixing of max-stable process based on their spectral representations. More recently, 
Kabluchko and Schlather [19] extend these results and obtain necessary and sufficient 
conditions for both mixing and ergodicity of max-i.d. random process. They define the 
dependence function of a stationary max-i.d. random process rj = {i](t)) t ez by 

„ , , Pfn(0) < a,ri(h) < a] „ . . , ^ 

rJh) = log ^'Y ~ ,L „\~ T i a > essinf n(0), h e Z. 
v ; 6 P[r/(0) < a]F[r](h) < a] /v 1 

Then, it holds with £ = essinf r?(0): 

• rj is mixing if and only if for all a > £, r a (n) — > as n — > +oo; 

• r] is ergodic if and only if for all a > £, n _1 Y^h=i T a(h) — > as n — >• +oo. 



Ergodicity is strongly connected to the strong law of large numbers via the ergodic theo- 
rem. The above results find natural applications in statistics to obtain strong consistency 
of several natural estimators based on non-independent but ergodic observations. 
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Going a step further, we address in this paper the issue of estimating the strong mixing 
coefficients of max-i.d. random fields. In some sense, ergodicity and mixing state that 
the restrictions rjs l and r]s 2 to two subsets S\,S2 become almost independent when the 
distance between Si and S2 goes to infinity. Strong mixing coefficients make this statement 
quantitative: we introduce two standard mixing coefficients a(Si, S2) and /3(Si, S2) that 
measure how much r]s 1 and i]s 2 differ from independence. The rate of decay of those 
coefficients as the distance between Si and S2 goes to infinity is a crucial point for the 
central limit theorem (see Appendix IA.4I) . As an application, we consider the asymptotic 
normality of three simple estimators of the extremal coefficients of a simple max-stable 
stationary random field on 7L d . 

Our approach differ from those of Stoev [25] based on spectral representations and of 
Schlather and Kabluchko [19j based on exponent measures. It relies on the Poisson point 
process representation of max-i.d. random fields offered by Hahn, Gine and Vatan [16] 
(see also Appendix IA. 1[) and on the notions of extremal and subextremal points recently 



introduced by the authors [14j. Palm theory for Poisson point process and Slyvniak's 
formula are also a key tool (see Appendix IA.3j) . 

The structure of the paper is the following: the framework and results are detailed in 
the next Section; Section 3 is devoted to the proofs and an Appendix gathers some more 
technical details. 



2 Framework and results 



Let (Q, J 7 , P) be a probability space and T be a locally compact metric space. We note 
C(T) = C(T, [0, +00)) the space of nonnegative functions endowed with the topology of 
uniform convergence on compact sets and C its Borel a-field. Let /1 be a locally finite 
Borel measure on Cq(T) = C(T) \ {0} satisfying 

£*[{/ G Cq(T); sup / > e}] < 00 for all compact K C T and e > 0, (1) 

and $ a Poisson point process (P.P.P.) on Co(T) with intensity \x. More rigorously, we 
should consider $ as a random point measure rather than as a random set of points, since 
there may be points with multiplicities. It is however standard to consider $ as a random 
set of points with possible repetitions. 

We consider the random process 

rj(t) = max{0(t), G $}, teT, 
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with the convention that the maximum of the empty set is equal to 0. Condition (CQ) 
ensures that the random process rj is continuous on T (see [16] and Appendix IA.ll) . 
Another property is worth noting: 77 is max-infinitely divisible. This means that for 
all n > 1, there exist independent and identically distributed continuous random fields 

{r]i,n)i<i<n such that 

where V stands for pointwise maximum and = for equality in distribution. Note that 
the max-infinite divisibility of r] is a simple consequence of the superposition Theorem 
for Poisson point processes. Furthermore, for all t G T, the essential infimum of the 
random variable r](t) is equal to 0. As shown by Gine, Hahn and Vatan |16| . up to simple 
transformations, essentially all max-i.d. continuous random process on T can be obtained 
in this way (see Appendix lA.il) . The measure n is called the exponent measure associated 



to the max-i.d. process r\. 

We now introduce the so-called a- and /3-mixing coefficients. For more details on 
strong mixing conditions, the reader should refer to the recent survey by C.Bradley [3] 
or to the monographs [T5| 1221 HJ |5j EJ [13]. For ScTa closed subset, we denote by J-$ 
the cr-field generated by the random variables {f]{s), s G S}. Let Si, S2 C T be disjoint 
closed subsets. The a-mixing coefficient (or strong mixing coefficient by Rosenblatt [23J) 
between the a-fields J r $ 1 and J-g 2 is defined by 

a(Si, S 2 ) = sup {\F(A HB)- F(A)F(B)\; A G J : s 1 , B G Ts 2 

The /3-mixing coefficient (or absolute regularity coefficient, see Volkonskii and Rozanov 
|27] ) between the cr-fields J r $ 1 and J"5 2 is given by 

P(Si,Sz) = \\P Sl US 2 -'Ps 1 ®Vs 2 \\var (2) 

= sup { \V Sl us 2 (C) - V Sl ® V s , (C)|; C G C Sl us 2 } 

where || • \\ var denotes the total variation of a signed measure and Vs is the distribution 
of the restriction r]s in the set C(5*) of continuous functions on S endowed with its Borel 
cr-field Cs- Since Si and S 2 are disjoint closed subsets, C(Si U S 2 ) is canonically identified 
with C(Si) x C(S 2 )- It is well known that 

1 1 J 

0(Si,S 2 ) = -supj^^inAn^O-P^P^Ol} 

i=i j=i 

where the supremum is taken over all partitions {Ai, . . . , Ai} and {Bi, . . . , Bj} of with 
the AiS in J 7 ^ and the B^s in J-5 2 . It also holds 

a(Si,S 2 )<^(3(Si,S 2 ). (3) 
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Our main result is the following. 

Theorem 2.1. Let r] be a max-i.d. process on T with exponent measure fi. Then, for all 
disjoint closed subsets Si, S 2 C T, 



P(S 1 ,S 2 )<2 [ nflt Sl V, f<s 2 v]^(df). 



In the particular case when Si and S 2 are finite or countable (which naturally arise for 
example if T = Z d ), we can provide an upper bound for the mixing coefficient /3(Si, S 2 ) 
involving only the 2-dimensional marginal distributions of the process rj. 

For (si,s 2 ) € T 2 , let /i SljS2 be the exponent measure of the max-i.d. random vector 
(rj(si), rj(s 2 )) defined on [0, +00) 2 by 

fi SUS2 (A) =»[{f e C (T); (f( Sl )J(s 2 )) eA}}, Ac [0,+oo) 2 Borel set. 
Corollary 2.1. If Si and S 2 are finite or countable disjoint closed subsets ofT, 

P(Si,S 2 ) < 2 I P M Sl ) ^ Vi> ^ ^ y2}^s 1 ,s 2 (dyidy 2 ). 

sieSi s 2 es 2 J 

Next, we focus on simple max-stable random fields, where the phrase simple means 
that the marginals are standardized to the standard unit Frechet distribution, 

F[r](t) <y] = expf-y- 1 ]!^}, yeR,teT. 

In this framework, an insight into the dependence structure is given by the extremal 
coefficients 9(S), S C T compact, defined by the relation 

P[su P 7?( S ) <y] = exp[-e(S)y- 1 ], y > 0. (4) 

ses 

Theorem 2.2. Let rj be a simple max-stable random field. 
For all compact S C T, 

C{S) = E[sup{^(s)- 1 ; seS}] < 00 (5) 

and furthermore: 

• For all disjoint compact subsets Si, S 2 C T, 

P(Si, S 2 ) < 2 \c(Si) + C(S 2 )] \e(Si) + 9(S 2 ) - 6(Si U S 2 ) 



5 



• For (Si,i)i e i and (S2j)jeJ countable families of compact subsets ofT such that S± = 
Uiej5'i ) i and S 2 = Uj £ jS 2 j are disjoint, 

In the particular case when Si and S*2 are finite or countable, the mixing coefficient 
j3(Si, S2) can be bounded from above in terms of the extremal coefficient function 

0{si,s 2 ) = 6({si,s 2 }), sx,s 2 eT. 

We recall the following basic properties: it always holds 9(si,S2) € [1>2]; 6(si,s 2 ) — 2 
iff Tj(si) and 7i(s 2 ) are independent; 6(81,82) = 1 iff t](si) = 77(52)- Thus the extremal 
coefficient function gives some insight into the 2-dimensional dependence structure of the 
max-stable field rj, although it does not characterize it completely. 

Corollary 2.2. Suppose 77 is a simple max-stable random fields. If Si and S 2 are finite 
or countable disjoint closed subsets of T, then 

P(S 1 ,S 2 )<4J2 Z>-0(ai,*O]- 

We illustrate our results on two classes of stationary max-stable random fields on R . 

Example 2.1. We consider the Brown-Resnick simple max-stable model (see Kabluchko 
et al. |18|). Let (Wj)j>i be independent copies of a sample continuous stationary in- 
crements Gaussian random field W = [W(t)) teR d with zero mean and variance <J 2 (t). 
Independently, let (Zi)i>\ be the nonincreasing enumeration of the points of a Poisson 
point process (0, +00) with intensity z~ 2 dz. The associated Brown-Resnick max-stable 
random field is defined by 

00 

7)(t) = \/ Ziexp[Wi(t) - a 2 (t)/2], t e R d . 

It is known that rj is a stationary simple max-stable random fields whose law depends 
only on the semi-definite negative function V, called the variogram of W, and defined by 

V(h) =E[(W(t + h) -W(t)) 2 ], heR d . 

In this case, the extremal coefficient function is given by 

6( Sl ,S 2 )=2V(y/V(s 2 -S 1 )/2), Si,S 2 CM*, 

where ^ denotes the cdf of the standard normal law. 



o(s 1)i ) + e(s 2 , j )-e(s 1>i \js 2ij ) . 
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Example 2.2. Our second class of example is the moving maximum process by de Haan 
and Pereira [12] , Let / : M. d — > [0, +oo) be a continuous density function such that 



/ f(x)dx — l and / sup fix + h) dx < 

JR d JR d \h\<l 



OO. 

h\< 

Let J^i=i $(Zi,Ui) be a Poisson random measure on (0, +oo) x W 1 with intensity z~ 2 dzdu. 
Then the random field 

oo 

V(t) = \/Z i f(t-Ud, tGR d , 
i=i 

is a stationary sample continuous simple max-stable random field. The corresponding 
extremal coefficient function is given by 

6(s 1 ,s 2 ) = / max(/(si — x), f(s 2 - x)) dx, si,s 2 eM. d . 
More generally, for any compact S C M d , 



= / max f( s ~ x ) dx. 

As noted in the introduction, our main motivation for considering the strong mixing 
properties of max-i.d. random fields is to obtain central limit Theorems (CLTs) for sta- 
tionary max-i.d. random fields. In this direction, we focus on stationary random fields 
on T = Z d and our analysis relies on Bolthausen's CLT [2J (see Appendix |A.4[) . 



We note \h\ = maxi^j^ \ hi\ the norm of h G Z d and \S\ the number of elements of 
a subset S C Z d . We note dS the set of elements h G S such that there is h' £ S with 
d(h, h') — 1. The distance between two subsets Si, S 2 C 7L d is given by 

d(Si, S 2 ) = min{|s 2 - si|; s x G Si, s 2 G 5*2}. 



A random field X = (X(t)) teZ d is said to be stationary if the law of (X(t + s)) t&I d 
does not depend on s G Z d . We say that a square integrable stationary random field X 
satisfies the CLT if the following two conditions are satisfied: 

i) the series a 2 = Yltez d Cov[X(0),X(t)] converges absolutely; 

ii) for all sequence A n of finite subsets of Z d , which increases to Z, d and such that 
lim^oo |<9A n |/|A n | = 0, the sequence |A n |~ 1/2 ^2 teAn (X(t) -E[X(t)]) converge in 
law to the normal distribution with mean and variance a 2 as n — > oo. 
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Please note we do not require the limit variance a 2 to be positive; the case a 2 = 
corresponds to a degenerated CLT where the limit distribution is the Dirac mass at zero. 
Bolthausen's CLT for stationary mixing random fields together with our estimates of 
mixing coefficients of max-i.d. random fields yield the following Theorem. 

Theorem 2.3. Suppose rj is a stationary max-i.d. random fields on Z, d with exponent 
measure \i and let 

7 (/i) = J Pfo(O) < yi, r){h) < y 2 ) ^ 0th {d yi dy 2 ), h e Z d . 

Let g : MP — » M. be a measurable function and t±, . . . , t p G Z d such that 
E[g( V (t l ),..., V (t p )) 2+5 ]<oo for some 5 > 0, 

and assume that 

oo 

l{h) = o^m 6 - 1 ) and ^ m* -1 sup 1 (h) s/{2+5) < oo. (6) 
\h\>m m=i \ h \^ m 

Then the stationary random field X defined by 

X{t)=g{r ] {t 1 + t),...,r ] {t p + t)), t G Z d 

satisfies the CLT. 

Condition Q requires that 7 goes fast enough to at infinity. It is met for example if 

l(h)<C-\h\- b for some b > dmax(2, (2 + S)/5) and C > 0. (7) 

If r\ is simple max-stable, we have ~f(h) < 2(2 — 6(0, h)), with 6 the extremal coefficient 
function. 

As an application, we consider the estimation of the extremal coefficient for a station- 
ary simple max-stable random field on Z d . For h G Z d , we note 6(h) = 6(0, h). Equation 
d3]) implies 

0(h) = -y]ogp(h,y) with p(h,y) = F(j)(0) < y, rj(h) < y), y > 0, 
suggesting the simple estimator 

8n\ h ) = -ylogp n (h,y) with p n (h, y) = |A n | _1 ^ l Wt)<y, r,(t+h)<y} (8) 

teAn 
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where A n is a sequence of finite subsets increasing to 7h d such that |<9A n |/|A n | — > as 
n — > oo. The fact that the naive estimator 9n\h) depends of the threshold level y > 
is not satisfactory. Alternatively, one may consider the following procedures. Smith pi] 
noticed that min(?7(0) _1 , 'q(h)^ 1 ) as an exponential distribution with mean 9{h)~ l and 
proposed the estimator 

0( 2 ) (h) = t^l 

Cooley et al. [7] introduced the F-madogram defined by 

u F (h)=E[\F( V (0))~F( V (h))\} with F(y) = exp(-l/y)l {y>0} 

and show that it satisfies 

/IN 19(h) -1 . , , „„, l + 2v F (/i) 

^ ) = 2^)Tl or eqmvalently W = __ 

This suggests the estimator 

n 1 ' \K\-2J2teA n \F(v(t))-F(r 1 (t + h))\- 

The following Proposition states the asymptotic normality of these estimators. 

Proposition 2.1. Suppose that r] is a stationary random field on 7* d with extremal coef- 
ficient function satisfying 

2 - 9(h) < C ■ \h\- b for some b > 2d and C > 0. (9) 

Then, the estimators 9 n (h), i — 1,2,3 are asymptotically normal: 

\K\ 1/2 {9 i n > (h)-9(h)) A/"(0, of) asn^oo 

u>zi/i lim^ variances 

a\ = y 2 J2( e Mmh)-9({0,h } t } t + h}))y- 1 ]-l), 

teZ d 

4 = ^(/^) 4 ^Cov[min(r ? (0)- 1 ,r / (/ i )- 1 ),min(r ? (t)- 1 ,r ? (t + / i )- 1 )], 

tez d 

a! = (W + l) 4 ^Cov[|F(r ? (0))-F(^))|,|F(r ? (t))-F(r / (t + / i ))|]. 



Interestingly, the function y \— > af is strictly convex, has limit +oo as y — >■ + or +oo 
and hence it admits a unique minimizer y* corresponding to an asymptotically optimal 
truncation level for the estimator §n \ Unfortunately, the limit variances o~\ and a\ are 
not so explicit so that a comparison between the three is difficult. 
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3 Proofs 



3.1 Strong mixing properties of extremal point processes 

In the sequel, we shall write shortly Co = Co(T). We introduce here the notion of S- 
extremal points that will play a key role in this work. We use the following notations: if 
fi, f 2 are two functions defined (at least) on S, we note 

fi =s h if and onl y if Vs G S, fi(s) = f 2 (s), 
fx < s f 2 if and only if Vs G S, < f 2 (s), 

fi i-s h if and onl y if 3s G S, fi(s) > f 2 (s). 

A point cf) G $ is said to be S'-subextremal if <fi <s tj, it is said S-extremal otherwise, 
i.e. if there exists s G S such that 4>(s) = i](s). In words, a S'-subextremal point has no 
contribution to the maximum r] on S. 

Definition 3.1. Define the S-extremal random point process $g and the S-subextremal 
random point process $^ by 

$<t = {(f) G $; 0^5^} and $ 5 = {<p G $; < 5 r?}. 

The fact that $jt and $^ are well defined point processes, i.e. that they satisfy some mea- 
surability properties, is proved in [H] Appendix A3. Clearly, the restriction rjs depends 
on <3>^t only: 

rj(s) = max{0(s); G $5}, s G S. 

This implies that the strong mixing coefficient (3(S\,S 2 ) defined by Equation (T5]) can be 
upper bounded by a similar /3-mixing coefficient defined on the level of the extremal point 
process $^ , $<t 9 . For % = 1, 2, let P $ + the distribution of $g. on the space of locally finite 

point measures on Co and let P (q> + $ + \ be the joint distribution of ($t , We define 

= ) - p ^ ® 11— ( io ) 

It holds 

/3(S 1; S 2 ) </?($+, $+). (11) 

The following Theorem provides a simple estimate for the /3-mixing coefficient on the point 
process level. It implies Theorem 12.11 straightforwardly and has a clearer interpretation. 

Theorem 3.1. • The following upper bound holds true: 

/?($+,$+) <2P[$+n$+=0] (12) 
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with 

H ®t = 0] < / P(/ V, f <s 2 17) Mf). (13) 

J Co 

• If the point process $ is simple (in particular in the max- stable case), the following 
lower bound holds true: 

W 1 ,$S,)>P[^n$S, = 0] (14) 

Clearly, equations (ITT]) . (I12p and (fT3"j) together imply Theorem 12.11 

Remark 3.1. The upper and lower bound in Theorem 13.11 are of the same order, and 
hence relatively sharp. It is not clear however how to lower bound f}(S\,S<i) and how 
sharp the upper bound in Theorem 12.11 is. 

3.2 Proof of Theorem I3TT1 

The upper bound for the mixing coefficient /3(<&g ,<&g ) defined by Equation (TlO|) relies 
on a standard coupling argument. There are indeed deep relationships between /3-mixing 
and optimal couplings (see e.g. chapter 5) and we will use the following well known 
result: 

Lemma 3.1. On a probability space (fl, J 7 , W), suppose that the random variables 
($^,$^) ; i — 1,2 are such that: 

i) the distribution of(^g^,^g^) is $+ ),' 

ii) the distribution of $| 2 2 ) is P$+ <g) P$+ . 

Then, 

/?($+,$+) <P 



We say that the random variables ($g , <&g ), i = 1,2 satisfying i) and ii) realize a 
coupling between the distributions Pt$+ $+ \ and P $ + <8> -P<j,+ • 

v Si' S 2 ' Si S 2 

In order to construct a suitable coupling, we need the following Lemma describing the 
dependence between and 

Lemma 3.2. Let S C T be a closed set. The conditional distribution of <3>g with respect 
to <3>g is equal to the distribution of a Poisson point process with intensity l{f <sr) y/j,(df). 
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Note that in the particular case when T is compact and $^ is finite almost surely, 
Lemma l3~2l follows from [14] Theorem 2.1 and Corollary 2.1. For T non compact, the proof 
need to be modified in a non straightforward way and is postponed to the Appendix IA.21 

We now construct the coupling providing the upper bound (fT2~l) . 

Proposition 3.1. Let (<&,rf) be an independent copy of (<&,rj) and define 

$ = $+ U {4> G $; 4> < Sl r,}. (15) 

The following properties hold true: 

• $ has the same distribution as $ and $ and satisfies 

K = K . = {0 e < a (16) 

• (^^j^sjj (^Sij^s 2 ) ^ s a coupling between P^+ $ + j and P^+ ®Pj>+ snc/i £/iat 



P 



Proof of Proposition \3.1\ 



^ <2P[<&+n<&£ = 0]. (17) 



Equation fl!6j) follows from the construction of $: consider 

Jftt) = \/ t G T; 

the maximum 77 is achieved on Si by the Si-extremal points $g , and the definition 
(fT5|) ensures that 77 and 77 are equal on S\ so that equation f fT6|) holds. 



Furthermore, conditionally on $J , the distribution of {0 G $; <5 a 77} is equal to 
the distribution of a Poisson point process with intensity 1{/< S v yfi(df). According to 
Lemma [3. 2 [ this is the conditional distribution of given $^ , whence ($5 , $5^) 
has the same distribution as ,$5 ). We deduce that $ = $^ U $g and $ = 
$5 U $5 have the same distribution. 



The coupling property is easily proved: since $ and $ have the same distribution, 



the law of (3?^, $s 2 ) is equal to P($+ )! since $ and $ are independent, ($5^ $s 2 , 



has law P$+ <g> P $ + . 



We are left to prove Equation f[T7|) . Since $J = $J , we need to upper bound the 
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probability P^^, 7^ ^sj- By construction, $ is obtained from $ by removing the 
points G $ such that -^Si V an d adding the points G 3>g . Hence, it holds 

{$+ ^ $+} C {30 G * a rf) U {30 G $+ , 77} 

Noting the equality of events 

{30 g $+, ^ v} = {30 e n $£} = {$£ n ^ 0}, 

we obtain 

Pft ^ *JJ < P[6+ n $+ ^ 0] + P[30 G rf 5l 77]. 
Since $ and $ have the same law, we have 

P[$+ n$+ ^ 0] = p[$+ n <&£ ^0]. 

Hence, equation ( TTTj) follows from the upper bound 

p[30g$+, 0^ 5l 77] < p[*+ n ^ 0] 

that we prove now. Using symmetry and exchanging the roles of $ and <3> on the 
one hand and the roles of Si and £2 on the other hand, it is equivalent to prove that 

p[30g$+, 0^<P[$+n$+ ^0]. 

We conclude the proof by noticing that the inclusion of events 
{30 G V} C {30 G 

entails 

P[30 g $+, <£s 2 77] < P[30 e r7] = n $£ 2 ^ 0]. 

□ 

We now complete the proof of Theorem 13.11 by proving Equations ( TT3]) and ( TT4l) . 
Proof of Equation f fTB"]) : We observe that 

{$+ n $+ ^ 0} = {30 g A, ^ 17} 

which entails 

p[$+n$+ ^0] <e[$^i 
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Noting that <f> -fts. 77 if and only if 
Appendix ) and compute 



-^Si max(f — {<p}), we apply Slyvniak's formula (see 



E [y^ max($- 



{4>}), 0^ S2 max(<f-{0})} 



E [l{/* Sl max(*), /^< ?2 max($)}M C (/') 



Co 



5^Si ^, / ^s 2 ^MO- 



□ 



Proof of Equation (fT4"j) : Clearly, for all measurable subset C C M p (Co) x M p (Co), 

> *£) e C] - P[(*+ , £+ ) g C]| 

where $ is an independent copy of We obtain the lower bound (JHJ) by choosing the 
subset 

C = {(Mi, M 2 ) G M P (C ) x M p (Co); M^M^ 0}. 

This yields indeed 

/?($+,$+) > |p[*+ n $£ 2 ^ 0] - P[*£ n $t 2 ? 0]|. 

and in the case when $ is a simple point process, i.e. when the intensity measure \i has 
no atom, we have 

p[$+ n ^ 0] = P[$ n $ ^ 0] = 0. 

□ 



3.3 Proof of Corollaries D and |2H and Theorem 12^1 

Proof of Corollary \2. li We have for all / G Co, 

{/ < Sl V, f <s 2 V} = {3(si, s 2 ) 6S1 x S 2 , > 7^1), /(s 2 ) > r](s 2 )} 

= U SieSl U S2e52 {?7(si) < f(si), rj(s 2 ) < f(s 2 )} 

whence, for Si and S 2 finite or countable, 

nf *a V, flts 2 v}<J2J2 ^(Si) < f(si), V (s 2 ) < f(s 2 )}. 

S1ES1 s 2 eS 2 
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As a consequence, the integral in Theorem 12.11 satisfies 

V, f £s 2 v] (*(df) 
< EE / Pfo(«i)</(«i), V(s2)<f(s 2 )] f i(df) 

V V / P[»?(si) < 77(^2) < y^Hs^dyxdyz). 



Co 



siSSi S2&S2 



In the last line, we have used the fact that fi SliS2 is the image of the measure \x under the 
mapping / ^ (/(si), /(s 2 ))- □ 



Proof of Theorem \2.2c We recall that for a simple max-stable random field, the exponent 
measure \x is homogeneous of order —1, i.e. fi(cA) = c~ l [i(A) for all A C Co Borel set 
and c > 0. Also the assumption of standard unit Frechet marginals implies 

flt(y) = y-\ teT, y > 0. 

These conditions implies (See Gine and al. [H] Propostition 3.2 or de Haan and Fereira 
J9] Theorem 9.4.1 and Corollary 9.4.2) that fi can be written as 

/»oo P 

K A ) = / 1 {rfeA} r' 2 dra{df) 

J0 JCo 

where a is a probability measure on Co such that 

fit) a(df) = 1 for all t G T, 



'Co 
and 

/ sup f(s) o-(df) < 00 for all compacts' C T. 
J Co seS 

Using this, note that for all compact S C T and y > 0, 

P[sup/(s) <y] = exp[-/i({/ e C ; sup f(s)>y})] 



exp 
exp 



l{su Ps£S rf(s)> y }r dra(df) 

Co 



y 1 / sup f(s) a(df) 

JCn S&S 



'Co seS 

It follows that the extremal coefficient 6(S) defined by is equal to 



9(S) = [ sup f(s)a(df). 

JCn s&S 
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We now consider equation (jHJ) and observe that for S C T compact, 



C(S) = E (inf 17(a)) 

seS 



-1 



so that we need to provide a lower bound for inf^s rj(s). To this aim, we remark that 

inf ri(s) = inf max0(s) > maxinf d>(s). 

ses ses 06* <t>e$ ses 

The right hand side is a random variable with unit Frechet distribution since 

P[maxinf <p(s) < y] = P[V0 G $, inf (p(s) < y] 
4>e$ ses 1 ses 



exp(-//({/ e C ; inf f(s) > y})) 



ses 



and 



/i({/eC ; mU(s)>y}) = y- 1 [ Mf(s)a(df). 

ses J Cq ses 

Hence, if J c inf s6 s f(s) a(df) > 0, we obtain 

(maxinf (b(s))' 1 ] = ( [ inf f(s) a(df )] < oo. 
s es YK JJ J \J Co seS JKJ W V 



C7(S)=E (inf tKs)) 



< E 



For arbitrary compact S C T, we may however have J c inf se s f(s) cr(df) = 0. But, 
if S = B(sq,e) is a closed ball with center sq and radius £, the monotone convergence 
Theorem implies 

/ inf f(s)a(df)^ [ f(s )a(df) = l as 0, 

so that J Co inf seB ( S0 £0 ) /(s) cr(df ) > and C(-B(s 0) £o)) < 00 f° r £ o small enough. The 
result for general S follows by a compacity argument: there exist Si, . . . , S/. and £i, . . . , £& 
such that 5 1 C U*_ 1 5(s i , £j). Hence, 



sup 77(5) 1 < max sup 77(5) 1 < Vj sup r/(s^ 



-1 



se5 



and 



C(5)=E[supr / ( S )- 1 l < Vfif sup ^H=VC(% £l ))< 
L «eS J ~l L seB(s . >£ .) J ' 



00. 



i=l 



This proves equation 
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The upper bound for /3(Si, S 2 ) given by Theorem 12.11 can be expressed as 



JCo 

= 2 [ [ P[r/ jt Sl V , rf ^ S2 r,]r~ 2 dra(df) 

JO JCo 

roo 

2 1 I E lu. 



JCo 

max 



2 / E 

'Co 



{r>inf SieSl ?7(si)//0i), r>inf S2e s 2 rj{s 2 ) / f{s 2 )} 

( inf ?H inf 7r4)i^/)- 



r 2 dra(df) 



neSi /(si)' ^ 2 e5 2 /(s 2 ) 
We then introduce the upper bound 



(19) 



max ( inf 



Si /(Si)' s% 2 f(s 2 )) 



< max ( sup T](si) , sup i](s 2 ) r ) min( sup /(si), sup f(s 2 )) 



whence we deduce 

/9(Si,&) 
= 2E 



max ( sup r)(si) 1 , sup 77(^2) x ) 

siG5i 



min(sup /(si), sup f(s 2 ))a(df) 

sieSi s 2 (es 2 



< 2E 



sup 77(51) 1 + sup 77(52) 1 



L S16S1 



s 2 es 2 



min(sup /(si), sup f(s 2 ))a(df) 

Co sieSi s 2 es 2 



= 2[C(^) + C(^ 2 )] [9(S 1 ) + 9(S 2 )-9(S l US 2 )}. 

In the last equality, we use Equation (jSJ) defining C(S') together with the following 
simple observation: in view of Equation (Tl8|) , the equality 

min(sup f(s x ), sup /(s 2 )) + max ( sup f(st), sup /(s 2 )) = sup f(s 1 ) + sup /(s 2 ) 
sieSi s 2 es , 2 sies 1 ! s 2 es , 2 sieSi S2e5 2 

yields after integration with respect to a(df) 



min( sup /( Sl ), sup f(s 2 )) a(df) + 9(S, U S 2 ) = 9(3,) + 9(S 2 ) 

Co sieSi S2 eS 2 



This proves the first point of Theorem | 
The second point is straightforward since 

tsi v, f <s 2 v] = Pp« el, f it Sl>i v, 3j e J, f A 2J . v\ 
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so that 

iei jeJ 

< 2^^[c(5' M ) + c(5 2 , i )] [e{s l>i ) + e{s 2>j )-e(s lji us 2>j )]. 

iei jeJ 

□ 

Proof of Corollary \2.2\ This is a straightforward consequence of the second point of 
Theorem 12.21 with S\ = U^^lsi} and 5*2 = U S2e s 2 {s2}- It holds indeed 

(3(s 1 ,s 2 )<2j2 E i c (M) + c (M)] WW) + ^})-%}uW)] 

fieSi S2&S2 

with 

6({ Sl }) = 6({s 2 }) = I , 6({ Sl }U{s 2 }) = e( Sl ,s 2 ) 

and 

C({ Sl })=C({s 2 }) = l. 

The last equality follows from the fact that, for all sg5, r](s) has a standard unit Frechet 
distribution and hence ^(s) -1 has an exponential distribution with mean 1. Hence we 
obtain 

0(Si,S 2 )<4X; £(2-0( Sl ,s 2 )) 

SlG^l S2GS2 

□ 



3.4 Proof of Theorems 12.31 and Proposition 12.11 

Proof of Theorem \2.3[ According to Bolthausen's CLT for stationary mixing random fields 
(see Appendix IA.4j) . it is enough to prove that the mixing coefficients akj(m) defined by 
Equation d2TJ) with X(t) = girjfa+t), . . . ,r}(t p +t)) satisfy Equations fl22|) . (ggj) and flMj) . 

For 5 C Z d , we define S 1 = Uf =1 {s + tj, sGS}. The inclusion of cr-fields 

a({X(s), seS})c a({ V (s), s G S}) 
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entails a comparison of the related a-mixing coefficients: for disjoint Si, 5*2 C Z, 

a x (S 1 ,S 2 )<a r '(S 1 ,S 2 ), 

where the superscript X or rj denotes that we are computing the a-mixing coefficient of 
the random field X oy 1] respectively. Furthermore, 

I Si | < p\Si\, i — 1,2 and d(S\, S 2 ) > d(Si, S 2 ) — A, 

with A = maxi<i<j<p d(ti, tj) the diameter of {t\, . . . , t p }. Hence, with obvious notations, 

<zH < a v pk Jm -A), MeNU {00}, m > 1. 

Corollary 12.21 implies 

o£,(m) < - |/i|) < p 2 fcZ sup 7 (t), fc, Z G N, m > 1, 

|t|>m-A 

and 

a x 00 (m)<al 00 (m-A)<p £ T (t), m > 1. 

|t|>m-A 

In view of this, Assumption ([6]) entails Equations ([22]) . (1231) and ( !24|) . so that the random 
field X satisfies Bolthausen's CLT. □ 

Proof of Proposition \2.1[ Let h G 7L d . We apply Theorem 12.31 to the stationary random 
field 

X(t) = l{ v (t)<y, r)(t+h)<y}, t G Z d . 

Clearly E[X(t)} = exp[-6(h)/y] and E[|X| 2+5 ] < 00 for all 5 > 0. Assumption © together 
with 7(t) < 4(2 — 6(t)) ensure that Equation ([6]) is satisfied for 5 large enough. Hence the 
estimator p n (h,y) is asymptotically normal: 

\K\ 1/2 (p n (h, y) - p(h, y)) =► A/-(0, ft 2 ) 

with limit variance 

PKv) = ^Cov[x(o),x(t)] 

= ^(exp(e({0,h,t,t + h})/y)-exp(29(h))/y)) > 0. 

The ^-method entails the asymptotic normality of the estimator 0^(h) = —ylogp n (h,y): 

\K\ 1,2 {ei{h)-9{h)) ^Af(0,al) 
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with limit variance 

a\ = y 2 exp(29(h)/y)fi = 2/ 2 ]T ( ex P[( 2 W " M, * + h }))/v\ ~ l ) ■ 

t&L d 

The proof of the asymptotic normality of On and 6n is very similar and we give only the 
main lines. Using Theorem 12. 3[ we prove that 

q n = \A n \- 1 min(r/(t)- 1 ,r/(t+/i)- 1 ) 

te|A„| 

is an asymptotic normal estimator of 6{h)~ l : 

with limit variance 

Pi = Gov [minMO)- 1 , vih)' 1 ), wm(n(t)-\ V (t + hy 1 )} . 

t€Z d 

(2) 

The 5-method entails the asymptotic normality of 9n (h) = l/q n (h) with limit variance 
Similarly, 

v F ^h) = \K\~ l \Hv(t))-F( v (t + h))\ 

te\A n \ 

is an asymptotic normal estimator of vp(h) = E[|F(?7(0)) — F(i](h))\]: 

with limit variance 

(3 2 3 = Y Cov[|F(i7(0)) - F( v (h))\, \F( V (t)) - F( v (t + h))\]. 

The ^-method entails the asymptotic normality of 

em = \±p§r 

1 - 2z/ jF>n (/i) 

with limit variance 

a! = (6(h) + l) 4 /3 3 2 . 

□ 
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A Auxiliary results 



A.l Structure of max-i.d. random processes 

The structure of sample continuous random processes on a compact metric space was 
elucidated by Hahn, Gine and Vatan [16j. We extend here their results for T a locally 
compact metric space, in order to cover the standard cases T = Z d or T = R d . Such 
extensions have been considered for max-stable models on M. (see |9j Chapter 9.6) but we 
have found no reference in the max-i.d. case. 

Let rj be a continuous max-i.d. random process on C(T, M). Define its vertex function 
h : T —> [—00, +00) by 

h(t) = essinf r](t) = sup{x G E; F(rj(t) > x) = 1}. 

We will always assume that h is continuous. We can then suppose without loss of gen- 
erality that h = 0. Indeed, if h is continuous and finite, we may consider rj — h which 
is a continuous max-i.d. random field with zero vertex function; and if h is not finite 
everywhere, we may consider exp(^) — exp(h) which is max-i.d with zero vertex function. 

We note C(T) = C(T, [0, +00)) the space of nonnegative continuous function on T 
endowed with the topology of uniform convergence on compact sets and Co(T) = C(T) \ 
{0}. 

Theorem A.l. • Let r\ = (r](t))t e T be a continuous max-i.d. process on T with vertex 
function h = 0. There exists a unique locally finite Borel measure on Co satisfying 
condition ([T]) ; called the exponent measure off], such that 



Uf =1 {/ G C ; f^) >y t } =- logP nJU MU) < Vi} 



for all k > 1, ti, . . . , tfe G T k and y\, . . . , y^ > 0. 

Conversely, for any locally finite Borel measure on Co satisfying condition (JT]), there 
exists a continuous max-i.d. process 77 onT with vertex function h = and exponent 
measure fi. It can be constructed as follows: let $ be a Poisson point process on Co 
with intensity \i and define 

7](t) = max{0(t), (f) G $}, teT. 



Proof of Theorem \A.l\ Let (T n ) n >i be an increasing sequence of compact sets such that 
T = U n >\T n . We suppose furthermore that T n is included in the interior set of T n+ \. 
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The space C(T) = C(T, M + ) of nonnegative continuous functions on T endowed with the 
topology of uniform convergence on compact sets can be seen as the projective limit of 
the sequence of spaces C(T n , IR + ) endowed with the topology of uniform convergence. For 
m > n > 1, we define the natural projections 

vr n : C(T) ->• C(T n ) and 7r n , m : C(T m ) -> C(T n ). 

For each n > 1, the restriction TT n (r]) = r]T n is a continuous max-i.d. process on the 
compact space T n and according to [16], there exists a locally finite exponent measure \i n 
on C (T„) = C(T n ) \ {0} satisfying equation 

^[uti {/ e c ;/(ti) > yi }] = -logpfn^i < y*}" 

for all > 1, ti, . . . , tfc G and yi, . . . , > 0. Furthermore, for all £ > 
Mn[«5 n , £ ] < oo where S n , £ = {/ G C(T n ); sup / > ej. 

Let n > 1 and e > be fixed. For n > n , define the finite Radon measure by 

tt°' £ [A] = fi n [A n tt-^.J, A C C(T„) Borel set. 
Clearly, the following compatibility conditions holds true: for m > n > no, 

r,n ,e _ ~n ,e -1 
rn r^m n,m - 

Theorem 5.1.1 in [1] state the existence of projective limit of Radon measures, it implies 
the existence of a finite Radon measure fx n °' e on C(T) such that 

It is then easily checked that the measure \x on Cq(T) defined by 

fjt[A] = sup{/i no ' £ [v4]; n > 1, e > 0}, Ac C (T) Borel set 
is locally finite and enjoys the required properties. □ 



A. 2 Proof of Lemma 13.21 



First, please note that the event {/ <s n} depends only on the restriction rjs and is hence 
measurable with respect to the cx-field generated by $^t. 

In order to prove the Proposition, let A\,...,Ak C Co be disjoint compact sets and 
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rix, . . . , rife > 0. Let A = U^ =1 Ai and n = Yli=i n i- We compute the conditional probability 
with respect to $^ of the event 

{$5(^1) = ni ,...,^ s (A k ) =n k }. 

This event is equal to {$ 5 G £} with B = {N E M P (C ); iV(Ai) = m, . . . , N(A k ) = n k }. 
We remark that that it is realized if and only if there exists a n-uplet [<f>\, . . . , <f> n ) of atoms 
of <3> such that: 

- the atoms <f>x, . . . , <f> n are S-subextremal; 

- the point measure $ — X^jLi ^<t>j nas 110 S'-subextremal atom in A, i.e. it belongs to 

D = {Ne M P (C ); #3 (A) = 0}. 



Then the n-uplet (0i, . . . , n ) is unique up to permutation of the coordinates. The above 
observations entail that for all measurable C C M p {Cq), 



n-l 



i=l 



Slyvniak's formula entails 

, s G C, $ 5 G B] 



E 



L {*+ e C} 1 {$ i (A)=0} n , 



,A 1 



(20) 



Summing this relation over the different values of rii, . . - ,n k G N and the related sets 
B = {N G Mp(C ); iV(Ai) = m, . . . , iV(A fc ) = ra fc }, we obtain 



P[$+ g C] = E 
So we can rewrite Equation fl20|) as 

P[$+ G C, $3 G B] = e" 
where 

exp [-//({/ G A; / <5"})] 



1 {*+gC} 1 {*J(A)=0} 



exp[u({/G A; f < s rj})]K(rj s , B) 



K(vs,B) 



nl 



is the conditional probability of {$3 G £>} with respect to $J (note it depends on $3 only 
through the restriction 775). We recognize the distribution of a Poisson random measure 
with intensity l{f <sV }fi(df) and this proves Lemma [3.21 
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A. 3 Slyvniak's formula 



Palm Theory deals with conditional distribution for point processes. We recall here one 
of the most famous formula of Palm theory, known as Slyvniak's Theorem. This will be 
the main tool in our computations. For a general reference on Poisson point processes, 
Palm theory and their applications, the reader is invited to refer to the monograph [26] 
by Stoyan, Kendall and Mecke. 

Let Mp(Co) be the set of locally-finite point measures iV on Co endowed with the 
cr-algebra generated by the family of mappings {iV i-> N(A), A C Co Borel set}. 

Theorem A. 2 (Slyvniak's Formula). 

Let $ be a Poisson point process on Co with intensity measure fi. For all measurable 
function F : Cq x M P (C ) -> [0, +oo) ; 

k k-1 
o i=l j=l 



E 



/ E[F(/i, ...,/*, $)] H m {dfi, ...,dfk)- 



A. 4 A central limit Theorem for weakly dependent process 

Since the pioneer work of Ibragimov [IT], many versions of the central limit Theorem for 
weakly dependent processes have been developed under various strong mixing conditions. 
We present here a central limit Theorem for stationary mixing random fields due to 
Bolthausen [2j. Let (Xk) keZ d be a real valued stationary random field and recall the 
definition of the a-mixing coefficient (l2Tj) . If A C Z d , we note |A| the number of elements 
in A and <9A the set of elements k G A such that there is I ^ A with d(k, I) = 1. Let A n 
be a fixed increasing sequence of finite subsets of Z d , which increases to Z d and such that 
lim^oo |0An|/|An| = 0. Let S n = EheA n ( x h - E[Jf k ]). 

Bolthausen's central limit Theorem is based on the mixing coefficients 

aki{m) = sup jc^Si, S 2 )] \St\ = k, \S 2 \ = I, d(Sx, S 2 ) > m| (21) 

defined for m > 1 and k, I G N U {oo}. 

Theorem A. 3. Suppose that the following three conditions are satisfied: 

«ioo(m) = o{m- d ); (22) 
J2m=i »7i d-1 a fc i(m) < oo for all k>l, l>\ such that k + I < 4; (23) 
¥,[\X h \ 2+& ] < oo and YZ=i m d ~ l \otxi{m)\ 5 tt 2+ V < oo for some 5 > 0. (24) 
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Then the series a 2 = J2hez d Cov[X , Xh] converges absolutely and if furthermore a 2 > 0, 
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